RAL-TR-1999-017 
December 1999 



Six-quark decays of off-shell W + W~ pairs 
via parton-level QCD 
in electron-positron annihilations 

S. Moretti* 

Rutherford Appleton Laboratory, 
Chilton, Didcot, Oxon 0X11 OQX, UK. 

Abstract 

We study the decays into six quarks of off-shell pairs of bosons produced in electron- 
positron annihilations, through the 0{a 2 s ) order in the strong coupling constant. We 
give explicit helicity amplitude formulae. We present numerical results in the context of 
phenomenological analyses of relevance at LEP2 and future Linear Colliders: such as Myy 
determinations, New Physics and Higgs searches. 
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1 Introduction 



Pairs of W bosons have copiously been produced at LEP2 and studied in great detail by the 
four experimental Collaborations over the past few years. In fact, one of the main goals of 
such collider is the determination of M\y with a target accuracy of 50 MeV or less. One of the 
detection strategies adopted to measure the W ± mass is the kinematic reconstruction of the 
W ± resonances through the momenta of their decay products, in the fully hadronic channel: 
W + W~ —>■ jets. Although the event reconstruction is made hard in this decay mode by the large 
number of tracks in the detector and by the usual uncertainties related to measuring jet energies 
and directions, a task much less complicated in the case of semileptonic decays, W + W~~ — ► 2 
jets ^ ± plus missing energy (with I — e, r), and despite the existence of theoretical biases due 
to the relatively unknown 'Bose-Einstein correlations' |l| and 'colour- rearrangement' effects, 
from which the semi-leptonic decays are immune, hadronic decays of W + W~ pairs have been 
successfully exploited at LEP2. On the one hand, they yield the largest decay rate. On the other 
hand, the knowledge of all momenta in the final state helps tightening the mass resolution. 

The problem with the fully hadronic mode is twofold. Firstly, because two identical decays 
take place in the same event, one has the phenomenon of mis-pairing of jets. That is, even in 
the ideal case in which all hadronic tracks are correctly assigned to the parton from which they 
originate, one has to cope with the ambiguity that it is in practise impossible to uniquely assign 
any pair among the four reconstructed jets to the parent W ± on the sole basis of the event 
topology. Of all possible combinations of di-jet systems, only one is correct. Thus, an intrinsic 
background exists in W + W~ —>■ 4 jet events, in terms of simple combinatorics. Secondly, because 
of the large hadronic multiplicity, one also has the phenomenon of mis-assignment of tracks. This 
is induced by the procedure adopted in selecting jets. This is generally done by resorting to 
so-called jet clustering algorithms ||, wherein the number of tracks is reduced one at a time by 
combining the two most (in some sense) nearby ones. The joining procedure is stopped by means 
of a resolution parameter, y cut , and the final 'clusters' are called jets. Here, the ambiguity stems 
from the fact that a track assigned to a cluster, the latter eventually identified as the parton 
originating from one of the W s, might have actually been produced in the fragmentation of 
another parton coming from the second decay. 

In both cases, the consequence is a distortion of the 'line shape' of the W ± resonance, which 
needs to be accurately quantified if one wants to achieve the foreseen precision in the W ± mass 
measurement. In order to estimate these effects, one can resort to phenomeno logical Monte 
Carlo (MC) programs (e.g., HERWIG g, JETSET/PYTHIA [| and ARIADNE f§). These 
represent a valuable instrument in this respect, as they are able to describe the full event, 
from the initial hard scattering down to the hadron level. On the other hand, Matrix Element 
(ME) models are acknowledged to describe the large angle distributions of the QCD (and QED) 
radiation better than the former (see, e.g., J7], which are in fact superior in the small angle 
dynamics. Whereas the above MC programs have been in use for long time, so that their features 
need not being recalled here, it might be useful to review at this stage the progress made in 
ME calculations of e + e~ — > W + W~ — > jets, as many of the achievements in this field are very 
recent . 

For a start, it should be mentioned that the amplitude for e + e~ — > W + W~ — > qqq'q' is very 
trivial to derive, in fact, more of a textbook example (here and in the following, q, q' and q" 



1 



refer to (anti) quarks produced in the W ± 1 W+ and g splittings, respectively). It represents the 
lowest-order (LO) contribution to the W + W~ — > 4 jet hadronic signal. Higher-order (HO) effects 
are those involving gluon emission: for example, the real one (i.e., a tree-level contribution) in 
e + e~ — > W + W~ —>■ qqq'q'g and e + e~ — ► W + W~ — > qqq'q'gg events, which have been calculated 
in Refs. || and (JCTT], respectively. One- loop QCD corrections to e + e _ — ► W + W~ — > qqq'q' are 
also known to date JTTJ , and they have been combined with the LO emission of Ref . || into the 
complete 0(a s ) result |TT . 



It is purpose of this paper to study the reaction e + e _ — > W + W~ — > qqq'q'q"q" at tree-level, 
through the order 0(ot\ m a 2 s ). Earlier accounts of this process, with emphasis on it being a 
background to top-antitop production and decay at a future Linear Collider (LC), were given 



in Ref. ||12|| . The complete set of Feynman diagrams needed to perform such a calculation can 
be found in Fig. |l]. The plan of the paper is as follows. In the next Section we highlight the 
phenomenological impact at LEP2 and a LC of six-quark production via W + W~ decays. In 
Sect. |3] we describe our method of computation using tree-level perturbative QCD. In Sect. f| 
we present our results whereas Sect. [5] contains our conclusions. 



2 Phenomenology of six-quark decays of W + W pairs 

In our opinion, there are a few good reasons to tackle the calculation of the process e + e _ — » 
W+W- -> qqq'q'q"q". 

As for Mw determinations, one should note the following. 

1. The availability of all HO corrections to the leading four-parton decay of W+ pairs is 
essential, both at tree and loop level, for the following reasons. On the one hand, from the 
point of view of perturbative calculations, it is evident that the 4-jet rate is constituted 
not only of the fraction of 4-parton events in which all partons are resolved (i.e., their 
separation is above the cut-off y cu t), but also by the (4 + m) parton configurations (with 
m > 1) in which m partons remain unresolved. On the other hand, because of the way 
an experimental 4-jet sample of W + W~ events is normally selected (see, e.g., fl3f), also 
the (4 + m) parton contributions with m jets resolved are relevant, as they naturally 
enter the candidate experimental sample of hadronic decays of W + W~ pairs. In fact, 
in order to maximise the event rate of the signal, one usually requires to reconstruct at 
ieast four jets, all with separation above a minimum y cnt . Only eventually these jets are 
forced into exactly four, by merging together those which are 'closer'. In this respect, one 
subtlety should be noted concerning W + W~ — > 6 quark events. The diagrams shown in 
Fig. [l] implicitly assume that the flavour of the quarks produced in the gluon splitting, 
labelled by 7 and 8, is different from that of any of the fermions generated by the W+ 
decays, indicated by 3, 4, 5 and 6. In fact, in case one or two quarks (or antiquarks) are 
identical, the numbers of diagrams doubles or quadruples, respectively, the new graphs 
being obtained from the old ones by exchanging one or two identical fermion legs, in all 
possible ways (a minus sign factorises too, for any of these exchanges). This follows from 
the fact that the flavour of a quark and its origin (i.e., whether it comes from a W+ or a g) 
are indistinguishable in the experimental hadronic sample. From the point of view of W+ 
mass studies, it is evident that this is source of two types of systematic effects. Firstly, 
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the total rate of W + W~ — > 6 quark decays naively obtained by summing over five flavours 
in the g* — > q"q" splitting (the approach used in parton shower models) could be mis- 
estimated. Secondly, additional distortion effects to the line-shape of the W ± resonance 
could occur. In the spirit of Refs. |], U, the reader should not object to calling such effects 
'flavour-rearrangement' or 'Fermi-Dirac correlations'. 

We further focus on possible New Physics and Higgs boson searches. 

2. The very fact that one of the two W s has a chance to decay into four jets reproduces at 
LEP2 and a LC a dynamics similar to that analysed in several instances at LEP1, when 
studies of 7*, Z — > 4 jet decays were performed^, in view of the possible existence of New 



Physics contributions due light gluinos g of the MSSM |fL4[| . The evidence of such effects at 
LEPI was mainly searched for in the context of the measurement of the three fundamental 
colour factors of QCD: Ca, C f (the Casimir operators of the fundamental and adjoint 
representations of the gauge group SU(N C ), respectively) and T F (the normalisation of 
the generators of the fundamental representation). In fact, under the assumption that 
SU(N C = 3) is indeed the gauge group of QCD, with N c the number of colours, a 
measurement of these parameters (in particular of Tr = N F T F ) can be converted into a 
constraint on the number of coloured fermions active at the energy scale at which the decay 
takes place. For example, N F would be increased (by approximately 3) from its SM value 
at LEP energies (i.e., N F = 5) by the additional presence of gluinos produced via a g — > gg 



splitting [p~5[| . In ordinary QCD, one gets Ca = 3 and C F = 4/3 (also, T F = 1/2). Since 
the results of these analyses have not excluded the possible existence of Supersymmetric 
(SUSY) events in LEPI data in the form of very light gluinos, with mass below 1-1.5 GeV 
[0, |T(J, we consider whether contributions of e + e~ — > W + W~ — > qqq'q'gg events possibly 
entering the six-jet sample produced at LEP2 and a LC can at all be disentangled. (Notice 
that the above mass region has escaped also the LEPI limits imposed through the running 



of a s , at the three-loop perturbative level [17|, in hadronic decays of heavy particles WE 



(e.g., Z and r), as well as those obtained from other experiments (KTeV, NA48, E761) 
searching for decays of gluino bound states^, see Plfl .) The much reduced W + W~ cross 
section at LEP2 and a LC, as compared to the Z one at LEPI, clearly disfavours such a 
possibility. In contrast, it would be intriguing to consider a selection strategy similar to 



that advocated in Ref. [22j, based on the detection of a secondary vertex possibly due to 



a gluino decay with lifetime between, say, 10 12 (the typical scale of b and c quarks) and 



1 Not quite, it could well be argued, given that in W + W~ — » 6 quark decays one has to first isolate a subset of 
four jets which come from the same decay, out of the original six, and since such an operation is in principle 
affected by the same mis-assignment problems already described. In practise, we will show that is rather easy to 
select such a subject, the latter preserving the typical angular properties of gauge bosons and fermions exploited 
in the experimental fits (see discussion later on). 

2 Some possible, Supergravity (SUGRA) inspired decay modes are R° — > 7r + 7r~7 and R p — > S a n + , where 
R° is the so-called glueballino bound state R° = (gg) and the 'photino' 7 is the lightest SUSY particle (LSP), 
whereas R p = (uudg) and S° = (udsg), with S° — > R°A. The mass(lifetime) of the R° is set by the theory to 
be in the range 1.4 - 2.2 GeV(10~ 5 - 10~ 10 sec) whereas that of the R p is 1.6 - 3.1 GeV(2 ■ lCT 10 - 2 • lCT 11 
sec). For alternative decay modes, in the Gauge Mediated SUSY breaking scenario, which have recently been 
proposed and not yet exploited experimentally, see Ref. |2(J. 
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(the coverage of the LEP detectors) sec. Besides, the fact that a W ± cannot decay 
directly into b quarks (apart from Cabibbo-Kobayashi-Maskawa suppressed channels, that 
we neglect here) implies that at LEP2 and a LC the overwhelming background from 
ordinary QCD due to e + e~ — > W + W~ — > qqq'q'gg events, with ^ b, should effectively 
be removed by asking for just a single vertex tagging displaying a decay lenght comparable 
or longer than that induced by b quarks, d& ~ 0.3 mm. This way, the surviving six-jet 
sample would only be composed by e + e~ — > W + W~ — > qqq'q'bb and e + e~ — > W + W~ — > 
qqq'q'gg events. These have comparable production rates, which is a most welcome result. 
After all, the smallness of the SUSY effect with respect to ordinary QCD dynamics was 
really the limiting factor of the experimental analyses performed at LEP1, given that 
even the very large next-to-leading-order (NLO) corrections to the four-jet sample leave 
the shape of the distributions used to fit Ca, Cf and Tp practically unaltered [f23|. For 
all the above reasons, we believe it then important looking into this aspect of six-jet 
phenomenology. 

3. Finally, we consider the possibility that six-quark decays of W + W~ pairs with a gluon 
splitting into bb pairs, i.e., e + e~ — > W + W~ — > qqq'q'bb events, can act as a background 
to Higgs signals possibly produced at a LC. In the Standard Model (SM) one has the 
production mechanism e + e~ — > </>Z, whereas in the Minimal Supersymmetric Standard 
Model (MSSM) one can have two, e + e~ — > $Z and e + e _ — > §A, where $ = H,h represents 
any of the two scalar Higgs bosons of the SUSY theory and A is the pseudoscalar one. If 
one recalls that Z bosons decay into bottom quark pairs some 15% of the times and that 
the SM Higgs decay channel </> — > W + W~ is the one with largest branching ratio (BR) 
for M,p ~ 140 GeV, then it is evident that a very large part of the SM Higgs signal at 
a LC would appear through a six-jet signature, further considering that each W ± boson 
decays hadronically with a 70% BR. Similarly, in the MSSMQ, the H -> W + W~ (of the 
heavy scalar Higgs) decay rate can be very large on a big portion of the (Ma, tan/?) plane 
(especially for Higgs masses below the top-antitop decay threshold and at low tan/5 |f24j] ) 
and so is the case of the h —>■ W + W~ (of the light scalar Higgs) one, though for a more 
restricted area of the MSSM parameter space. Furthermore, A —>■ bb decays of the MSSM 
pseudoscalar Higgs are always dominant, except at low tan/5 and large Ma (above 200 
GeV). Indeed, considering that in the fully hadronic channel, in order to suppress both 
the effects of combinatorics and of reducible backgrounds (e.g., e + e~ — > 7*, Z* — > qqgggg 
|p5fl ), vertex tagging techniques will certainly be exploited to select bb pairs produced in Z 
and A decays, then it is quite likely that six-quark decays of pairs (with two 6's) can 
be a serious noise. Even more so if one further realises that these are irreducible in two 
respects: not only in the particle content of the final state, but also because they naturally 
contain a di-jet pair resonating at the mass. 



3 Assuming that the mass scale of SUSY partners of ordinary matter is in the TeV range, as preferred by 
many Grand Unification Theories (GUTs). 
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3 Calculation 



In the numerical part of our calculations, as centre-of-mass (CM) energies representative of 
LEP2 and a future LC, we have used the values E cm = 172 and 350, 500 GeV, respectively. 
As for the parameters of the theory, we have adopted Mz = 91.17 GeV, Tz = 2.516 GeV, 
Mw = 80.23 GeV, IV = 2.2 GeV, sin 2 9w = 0.23, a em = 1/128 and the two-loop expression 
for a s , with A^^ 5 = 0.200 GeV. Furthermore, we have kept all quarks massless as a default 
except the bottom ones (for which we used mj = 4.95 GeV), in order to speed up the numerical 
evaluations. (Electron and positron have mass zero too, so have the neutrinos.) In practise, 
as we neglect Cabibbo-Kobayashi-Maskawa mixing terms (see below, Tab. |3|), this corresponds 
to always neglect the mass along the fermion lines connected to the W ± currents in Fig. ITL 
while keeping the one of those emerging from the gluon. This should be a good approximation. 
On the one hand, the energy produced in the W ± decays both at LEP2 and a LC is typically 
much larger than the quark masses. On the other hand, the latter can give sizable effects in the 
splitting of the gluon [p6fl . This approach also naturally allows us to study effects of massive 
gluinos. As for these, we have spanned their mass m- g over the range to 10 GeV. While doing 
so, we also have to compare the six-quark and four-quark-two-gluino rates with the four-quark- 
two-gluon ones. That is, we have to calculate the process e + e~ — > W + W~ — > qqq'q'gg. To do 



so, we have resorted to the HELAS [£7j subroutines, whose results agree with those reported in 
Ref. [0. 

In addition, when e + e~ — > Higgs Z —> W + ^W~^bb — > bb + 4 jets and e + e~ — > $ A —> 
W + (*^W~^bb — > bb + 4 jets have been calculated (at tree-level), where Higgs = </>,$, with cf> 
the SM scalar boson and $ = H, h the corresponding ones in the MSSM, we have used again 
heavy b quarks, with their mass appropriately run up to the mass of the decaying pseudoscalar 
A boson in the second process, to match the procedure employed to compute the total widths 
of the latter J24|. Mass relations and couplings involving the neutral MSSM Higgs bosons 



have been computed using the full one-loop and the leading two-loop corrections |p8[| , using 
tan (3 and Ma as inputs, assuming a universal soft Supersymmetry-breaking mass of 1 TeV and 
negligible mixing in the stop and sbottom mass matrices. To calculate the Higgs cross sections 
and differential rates we have used the exact 2^6 MEs, including finite width effects of all 
unstable particles, that we have produced using again the HELAS subroutines. 

In order to calculate the Feynman diagrams in Fig. [I] we have used two different spinor 
methods, this enabling us to check the correctness of our results. The first one is based on the 



formalism of Ref. |p9|| . A second approach is based on the method developed in Refs. |30], [3T], [32 
Since we will express the helicity amplitudes in this last formalism, we will devote same space 
here to describe its technicalities. 

In this method, all spinors for any physical momentum are defined in terms of a basic spinor 
of an auxiliary light-like momentum. By decomposing the internal momenta in terms of the 
external ones, using Dirac algebra and rewriting the polarisation of external vectors by means 
of a spinor current, all amplitudes can eventually be reduced to an algebraic combination of 
spinors products u{pi)u(pj) (with ... labelling the external particles). 

(i) Spinors. Since all vector particles entering our calculation eventually splits into fermions, 
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AiA 3 


X{pu X 1 ;p 2 ;p 3 , ^3;c R ,c L ) 


++ 


(/il^2 + yU2^l) (Cfi/i 2 ^3 + c L n 3 r] 2 ) + c R S ( + , pi, p 2 )<S'(-, p 2 , P3) 


+- 


cl(hiV2 + A i 2^i)5'(+,P2,P3) + c L (c L fi 2 r]3 + c R fi 3 r]2)S(+,p 1 ,p 2 ) 



Table 1: The X functions for the two independent helicity combinations in terms of the functions 
5*, 1] and defined in the text. The remaining X functions can be obtained by flipping the sign of 
the helicities and exchanging + with — in the S functions and R with L in the chiral coefficients. 



A1A2A3A4 


Z(pi, Ax;p 2 , A 2 ; p 3 , A 3 ; p 4 , A 4 ; c R , c L ; c' R , c' L ) 


+ + ++ 


-2[S(+,p 3 ,pi)S(-,p 4 ,p 2 )c' R c R - ^i^V^c'rCl ~ Hic' L c R } 


+ + +- 


-2rj 2 c R [S(+,p 4 ,px)ii 3 c' L - 5 , (+,p 3 ,pi)/i 4 cy 


+ + -+ 


-2r) 1 c R [S(-,p 2 ,P3)fi4c' L - 5 , (-,p 2 ,p 4 )/i3cy 


+ - ++ 


-2r] 4 c' R [S(+,p 3 ,pi)fM 2 c R - 5 , (+,p 3 ,p 2 )/iiC L ] 


+ + — 


-2[5 , ( + ,Pi,p 4 )S , (-,p 2 ,p 3 )c' L C R - {ll/J>2V3V4C' L C L - VlV2^4C' R C R ] 


+ - +- 





+ - -+ 


-2[/i 1 /i 4 r] 2 ?7 3 c^c L + n 2 ^3Vmc' R c R - y u 2 /i 4 ^i^ 3 c / i c i? - /^/^r^c^Ci] 


+ 


-2rj 3 c' L [S{+,p 2 ,p 4 )^ic L - 5'(+,pi,p 4 )yU 2 c R ] 



Table 2: The Z functions for all independent helicity combinations in terms of the functions S, 
r) and fi defined in the text. The remaining Z functions can be obtained by flipping the sign of 
the helicities and exchanging + with — in the 5* functions and R with L in the chiral coefficients. 



we only need to introduce the treatment of the spinor fields. External fermionsQ of mass m and 
momentum p^ are described by spinors corresponding to states of definite helicity A, u(p, A)Q, 
verifying the Dirac equations 

^u(p, A) = ±mu(p, A), w(p, A)|/ = ±mu(p, A), (1) 

and the spin sum relation 

U (P> A MP' A) = p 7 ± m, (2) 

A=± 

where the sign +(— ) refers (here and in the following) to a particle (antiparticle). One can choose 
two arbitrary vectors k$ and k\ such that 

k ■ k = 0, ki ■ k\ — —1, k ■ ki — 0, (3) 

and express the spinors u(p, A) in terms of chiral ones w(k , A) as 

u(p, A) = w(p, A) + fxw(k , -A), (4) 

4 We shall use the term 'fermion' and the symbol 'u' for both particles and antiparticles. 
5 Here, p(X) represents a generic (anti)spinor four-momentum(helicity). 
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where 

w(p,X) = ^w(k ,-X)/r], (5) 



and 



The spinors w(ko, A) satisfy 



and therefore 



rn 



fi = ±—, r]=^/2\p-k \. (6) 

w(k , X)w(k , A) = 1 + ^ l5 ^ Q , (7) 



5^ w(ko, \)w(ko, A) = U. (8) 



The phase between chiral states is fixed by 

w(k ,\) = \k / w(k ,-\). (9) 

The freedom in choosing k$ and k\ provides a powerful tool for checking the correctness of any 
calculation. A convenient, though not unique choice, is the following: k Q = (1,0,0,-1) and 
k\ = (0, 1,0,0). In such a case the massless spinors in the two methods J29| and |30) coincide 
exactly, so that it is possible to compare in greater detail the two corresponding numerical codes. 
In particular, the results obtained with the two formalisms must agree for every single diagram 
and every polarisation of external particles. 

(ii) The S , X and Z functions. Using the above definitions one can compute the spinor 
functions 

S(X,p u p 2 ) = [u(pi, A)«(p2, -A)], (10) 
X(p 1 ,Xi;p2;P3,X 3 ;c r ,cl) = [u(pi, Ai)^4r«(p 3 , A 3 )], (11) 

Z(pi, Xi,p 2 , A 2 ; p 3 , A 3 ; Pa, A 4 ; c r , c l ; c' r , c' l ) = 

[u( Pu Aa)r"«(p2, A 2 )] [u(p 3 , A 3 )I>(p 4 , A 4 )], (12) 

r (')/i = 7 M r (') ) (13) 

rn = c<2p R + $p L , (14) 

the chiral projectors. 

By computing the resulting traces one easily finds (e 0123 = 1 is the Levi-Civita tensor) |3(J j3ll 



and 

where 

and 

with 



C / . x o (Pi ■ k o)(P2 ■ h) - (pi ■ fci)(p 2 ■ fco) + ^^poKKVxVl 

b{+,p 1 ,p 2 ) = 2 , (16) 

V1V2 
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for the S functions and the expressions listed in Tabs. [I] and |2] for the X and Z functions, 
respectively. For the S functions, one has 5'(— ,pi,p 2 ) — S(+,p 2 ,Pi)*, while the remaining X 
and Z functions can be obtained as described in the captions of Tabs. 

Other than the spinor parts, to each of the basic amplitudes are associated propagators 
functions. In the case of off-shell fermions, they have the form 

D ' ( V K) = en^W (17) 

whereas in the case of bosons one gets 

i [EiPi) - Mfr + iM v T v 

In eqs. (|T7|)-(|T8|), / is the flavour q, q' of a virtual fermion line, whereas V = W ,j,Z or g, 
being My = T v = if V = 7 or g. 

As for the couplings, the notation (c^ V , c~p V ) will refer to the pair of chiral indices (cr, c£) 
of eq. (0) entering in the expressions given in Tabs. and associated with the vertex involving 
a fermion / = e,q, q', q" (whose label will only appear if the vertex is flavour dependent) and 
a gauge vector V = W ± ,'j, Z, g, according to Tab. [| (Note that we will only need using 
the flavour- dependent terms (cfj", cf^), i.e., for / = e and V = J,Z.) For convenience, we also 
introduce the relative couplings g 1 — 1 and gz = g 1 j ta.n9 w entering the r yW + W~ and ZW + W~ 
vertices, respectively, where 6 W is the Weinberg angle. 
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z 


W ± 


9 


CR 


Q f 


9 R / swcw 





1 


c l 


Q f 


gi/s w c w 


l/V2s w 


1 



Table 3: The couplings cr and cl of eq. fll4f) for u and d type (anti)quarks and electrons/positrons 
to the gauge bosons 7, Z, W ± and g. One has (adopting the notations sw = sm9w and 
c w = cos9 w ) 9r = ~Q f s 2 w and g{ = Tl - Q f s 2 w (with q = u,d), where (Q U ,T£) = (+|,+|), 
(Q d ,T-^) = (— |, — |) and (Q e ,T 3 e ) = (—1, — |) are the fermion charges and isospins. 



We are now ready to present the explicit expressions of the helicity amplitudes associated 
to our process. To do so, we conventionally assume that initial state momenta are incoming, 
whereas the final state ones are outgoing. This way, we can define —hi — bj — 1, where i — 1, 2 
and j = 3, ...8, so that J2k=i ...s^kPk — 0. In correspondence to the graphs in Figs. [3], one can 
write the Feynman amplitude squared, summed/averaged over final/initial colours and spin, as 

\ M \ 2 = ^EEE^.({^(W), (19) 

{A} i=l j=l 

where Cy are the colour factors (see below) and with referring to a summation over all 
possible combinations of the helicities Ai, ...As of the external particles. The quantities g s and 
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e are related to the aforementioned couplings by the usual relations g\ = Ana s and e 2 = 47ra ( 
(in natural units). Assuming, for sake of illustration, that the process (via W + W~) 



e+(p!,Ai) e (p 2 ,A 2 ) — > u(p 3 , A 3 ) rf(p 4 , A 4 ) s(p 5 , A 5 ) c(p 6 , A 6 ) b(p 7 ,X 7 ) fr(p 8 ,A 8 ) 
has to be calculated, then the helicity amplitudes Tj can be written as 

^({A}) = -AXP7+P8)£W^3+P4 + P7+P8)IMP5+P6) (20) 
A*(P3 + P7 + Ps)D v& {P2 -Pb~ Pd) 

E E &AEE 

i=3,7,8 j=2, 5,6 A=±A'=± 

Z(p 7 , A 7 ; p 8 , -A 8 ; Ps, A 3 ; Pi, A; c^, c£; 4, cf ) 
^(Pi, A; p 4 , -A 4 ; Pi, Xi,Pj, A'; c^, cf ; c^, cf ) 
Z{Pj,\;P2,-M;P5,M;P6,-\;c R ,c L ;c R ,c L ), 



T 2 ({A}) = +D g (p 7 +p 8 )D w (p 3 +p 4 +p 7 +p 8 )D w (p 5 +p 6 ) (21) 

Ai(P4 + P7 + p 8 )D Ue (p 2 -p 5 - P&) 

E E &AEE 

i=4,7,8 j=2,5,6 A=±A'=± 

Z(p 7 , A 7 ; Ps, -A 8 ; Pi, A; p 4 , -A 4 ; 4, c£; cfj, c|) 
Z\Pz, A 3 ; ^, A; p 1: Xi,p j: X ;c R ,c L ;c R ,c L ) 
Z{Pj,X;p2,-X 2 ;p 5 ,X 5 ;p 6 ,-XQ;c R ,c L ;c R ,c L ), 

^3({A}) = +D g (p 7 +p 8 )D W (p 3 +p 4 )D W (p 5 +p 6 +p 7 +p 8 ) (22) 
£>s(P5 + P7 + Ps)D Ue (p 1 - p 3 - P4) 

E E &AEE 

i=5,7,8 j=l,3,4 A=±A'=± 

Z(p 7 , X 7 ; ps, -A 8 ; p 5 , A 5 ; p*, A; cfj, cf ; cfj, c£) 
Z(pi, A; p 6 , -A 6 ; p j? A'; p 2 , -A 2 ; c^, cf ; cf , cf ) 
Z(pi, Ai;pj, A'; p 3 , A 3 ; p 4 , -A 4 ; cf , cf ; cf , cf ), 

T 4 ({A}) = -A 7 (p7+P 8 )-Dw(P3+P4)-tV(P5+P6+P7+P8) ( 2 3) 
£>c(P6 + P7 + Ps)Du e (Pl - P3 - Pi) 

E E to* E E 

j=6,7,8 j=l,3,4 A=±A'=± 

^(P7, A 7 ; Ps, -A 8 ; Pi, A; p 6 , -A 6 ; cfj, c£; c^, c£) 
Z{p 5 , A 5 ; Pi, A; pj, A ; p 2 , -A 2 ; c fl , c L ; c fl , c L ) 
^(Pi, Ai; Pj, A'; p 3 , A 3 ; p 4 , -A 4 ; cf , cf ; cf , cf ), 
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^({A}) = +D g (p7+p 8 )D w (p 3 +p4+p 7 + p s )D w (p 5 +p 6 ) (24) 
D U (P3 +P7+ Pa) E 9vD v {pi + Pn) 

E biE 

i=3,7,8 A=± 

{Z(pu Ai;p 2 , -A 2 ; Pi, A; p 4 , -A 4 ; cf , cf ; cf , cf ) 
[ E & i x (P5,A 5 ;Pj;P6,-A 6 ;cR ,cf ) - 

.7=3,4,7,8 

E & i x (Ps, A 5 ;Pj;p 6 , -A 6 ; cf , cf )] + 

J"=l,2 

^(Pi, A; pa, -A 4 ; p 5 , A 5 ; p 6 , -A 6 ; c^ , c L ; c R , c L ) 

[ 53 h j X {PlAl]Pj]P2,-X 2 ]cf ,cf ) - 
3=5,6 

53 b j X(p 1 , X 1 ;p j ;p 2 , -A 2 ; cf , cf )] + 

.7=3,4,7,8 

^(Ps, A 5 ; p 6 , -A 6 ; Pi, A x ; p 2 , -A 2 ; cf , cf ; cf , cf) 
[ 53 bjX(pi, A; pj-, p 4 , -A 4 ; cf , cf ) - 

J"=l,2 

53 b j X(p i ,X;p j ;p 4 , -A 4 ; cf , cf )]} 

.7=5,6 

Z(p 7 , A 7 ; Pa, -A 8 ; p 3 , A 3 ; p i; A; cfj, c£; c^, c£), 



^6({A}) = -^(P 7 +P 8 )^(P3+P4+P7+P8)^(P5+P6) (25) 

Az(P4 + P7 +Ps) E 9vD v (p 1 +p 2 ) 

V=~/,Z 

E 6<E 

i=4,7,8 A=± 

{Z(pu Ai;p 2 , -A 2 ; p 3 , A 3 ; p*, A; cf , cf ; cf , cf ) 
[ E 6 j- x "(p5,A 5 ;p j -;P6,-A6;cb , cf ) - 

,7=3,4,7,8 

E & j'A:(p5, X 5 ;pj;p6, -A 6 ; cf , cf )] + 

3=1,2 

Z{P3, A 3 ; p^ A; p 5 , A 5 ; p 6 , -A 6 ; % , c L ; c R , c L ) 
[ 53 6jX(pi, X 1 ;p j ;p 2 , -A 2 ; cf, cf ) - 

.7=5,6 

53 6i-X"(pi,Ai;pj;P2,-A 2 ;c^,Ci V )] + 

.7=3,4,7,8 

Z(p 5 , A 5 ; p 6 , -A 6 ; Pi, A x ; p 2 , -A 2 ; cf , cf ; cf , cf ) 
[ E & iAT(p 3 , A 3 ; pj-,pi, A; cf , cf ) - 

3=1,2 
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E b j X (P3,^Pj]Pi,^,CR,cY)]} 
j=5,6 

Z(p 7 , A 7 ; Ps, -A 8 ; Pi, A; p 4 , -A 4 ; Cr, (f L ; 4, cf,), 



T 7 ({\}) = -D g {p 7 + p 8 )D w (p 3 +p 4 )D w (p 5 +p 6 +p 7 +p 8 ) (26) 

A.(P5 +P7+ Ps) E 9vDvi.Pi + p-i) 
v=i,z 

E fcE 

i=5,7,8 X=± 

{Z(p u Xi,p 2 , -A 2 ; p i? A; p 6 , -A 6 ; ; c%, cf ) 

[ E b j x (Ps,h;Pj;P4,-X 4 ;c^,cY ) - 

.7=5,6,7,8 

E A 3 ; Pj-; p 4 , -A 4 ; c R , cf )] + 

.7=1,2 

Z(pi, A; p 6 , -A 6 ; p 3 , A 3 ; p 4 , -A 4 ; c^, cf ; c^, cf) 

[ b 3 X (PU X UP3^P2,-^2;cf,cf)- 
.7=3,4 

b j X(p 1 ,X 1 ;p j ;p 2 ,-X 2 ;c R v ,cf)} + 

i=5,6,7,8 

^(P3, A 3 ; p 4 , -A 4 ; pi, A x ; p 2 , -A 2 ; c^, cf ; cf, cf) 
[ bjX(pi, A; pj; p 6 , -A 6 ; c^, cf ) - 

J"=l,2 

X] b j X(p i ,X;p j ;p 6 , -A 6 ;c^,cf)]} 

.7=3,4 

Z(p 7 , A 7 ; Ps, -A 8 ; p 5 , A 5 ; Pi, A; cfj, c£; c^, c£), 



^s({A}) = +D g (p 7 +p 8 )D w (p 3 +p 4 )D w {p 5 +p 6 +p 7 +p 8 ) (27) 
^ C (P6 + P7 + Ps) E 9vD v (pi + P2) 

V=~f,Z 

E feiE 

i=6,7,8 A=± 

{Z(pi, Ai;p 2 , -A 2 ; p 5 , A 5 ; p i? A; c^ y , ; c^, cf ) 

[ E & i^(P3,A 3 ;p i ;p 4 ,-A 4 ;c^,cf ) - 
.7=5,6,7,8 

E & i x (^3, A 3 ; Pj-; p 4 , -A 4 ; c^, cf )] + 
3=1,2 

Z[P5, A 5 ; Pi, A; p 3 , A 3 ; p 4 , -A 4 ; c fi , c L ; c R , c L ) 
[ b j x (Pi> Ai;Pj;p 2 , -A 2 ; cff, ) - 

J=3,4 
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Z bjXipx, Xi;pj;p 2 , -A 2 ; cf , cf)] + 

j=5,6,7,8 

Z(P3, A 3 ; p 4 , -A 4 ; Pi, \i,P2, -A 2 ; c™, ; cf , ) 
[£ & i x (^5, A 5 ;Pj;Pi, A;cf ,cf ) - 

Z & i x (p 5 , A 5 ;pj;pi, A;c« ,cf )]} 

.7=3,4 

Z(p 7 , A 7 ; p 8 , -A 8 ; A; p 6 , -A 6 ; cfj, 4; cfj, 4). 

Concerning the colour factors, there are basically only two of these, if the flavours in the 
final state are all different. They are (hereafter, Nc = 3): Cy = —F-(Nq — 1) = 6, if i and j are 
diagrams in which the gluon emission takes place from the same W^, and Cy = otherwise 
(because of colour conservation )[]. 

If one or two quark flavours in the final state are identical, then the number of Feyn- 
man graphs proliferates, as explained in Sect. []. The spinor part of the additional dia- 
grams can easily be obtained by interchanging the labels of identical particles in the above 
formulae, factorising a minus sign for each of these operations. In addition, as many fac- 
tors of the form ^ multiply the amplitude squared as the number of n-tuple of identical fi- 
nal state particles. As for the new colour factors, one has to do some more work. However, 
it is rather trivial to realise than only the following colour structures need to be computed^: 

(^^j =N C = 3 



-O-Q-Q-Q 
LQJLQJJ. 




6 In other terms, 'perturbative' colour-rearrangement is not possible in W + W~ — > 6 quarks at 0(a 2 s ), contrary 
to the case of four-quark- two-gluon decays at the same order |10|. 

7 The product of the first two yields the factor dj = 6 mentioned above whereas the third produces the other 



Ci 



0. 
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= 



Since the all procedure is quite cumbersome, we refrain here from building up explicitly the 
correct MEs for identical flavours by combining the above colour factors (and their products) 
with the appropriate interferences among spinor amplitudes. Instead, we make available upon 
request our programs, that do include the described implementation. 

Before proceeding further, we would now like to devote some space to describe the procedure 
adopted to integrate the squared amplitude in eq. (p!9[). In fact, in order to deal numerically with 
the non-trivial resonant structure of our six-quark process, one has to apply some special care. 
Here, we have adopted the technique of splitting the ME in a sum of non-gauge-invariant pieces, 
each of these implementing a different resonant structure, and of integrating them separately 
with the appropriate mapping of the phase space variables. 

Things go as follows. Firstly, one isolates the diagrams with similar resonant structure by 
grouping these together in 'subamplitudes'. From the graphs in Fig. [TJ, one can recognise the fol- 
lowing two resonant structures: say, (a) W + — > (3478) and W~ — > (56) (graphs 1,2,5,6, so that 
T a = E^i.2,5,6 (b) W + -> (34) and W -> (5678) (graphs 3,4,7,8, so that T b = £ i=3i 4,7,8 
Secondly, one defines the mentioned non-gauge-invariant components of the amplitude squared, 
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by appropriately combining the subamplitudes. For example, we have simply taken the square 
of the two resonant subamplitudes and their interference: |T a | 2 , |T b | 2 and 2 Real (T a T b *), re- 
spectively Thirdly, one maps the phase space around the resonances. Fourthly, the various 
amplitude squared terms are integrated separately and added up in the end (to recover gauge 
invariance) to produce total and differential cross sections. In this respect, we would like to 
mention that all results presented here have been obtained by resorting to the adaptive multi- 
dimensional integrator VEGAS [33], and they have been counter-checked against the outputs of 



the multi-particle phase space generator RAMBO |34 



4 Results 

In order to select a six- 'jet' sample we apply a jet clustering algorithm directly to the 'quarks' 
in the final state of e + e~ — > W + W~ — > qqq'q'q"q" ' . For illustrative purposes, we use the Durham 
jet-finder |35j only. However, we remark that none of the main features of our analysis depends 



drastically on such a choice. This algorithm is based on the '(squared) transverse- momentum' 
measure 

2min(E 2 ,E 2 )(l-cos%) 
Vij = ~ , (28) 

where Ei and Ej are the energies and 9ij the separation of any pair ij of particles in the final 
state, with 3 < i < j = 4, ...8, to be compared against a resolution parameter denoted by y cut . 
In our tree-level studies, the selected rate is then nothing else than the total partonic cross 
section with a cut j/^ > y cut on any possible ij combination. 

Fig. |2] presents the total cross section at LEP2 (above) and a LC (below) for six-quark 
events as a function of the resolution, with and without the correlations described in Sect. ^. 
(A summation over all possible combinations of quark flavours has been performed.) We see 
a large effect on the integrated rates. Indeed, the ratio between the two curves is about 1.6 
for y cnt ~ 0.01 at both energies. The cross section including the interference effects is indeed 
smaller, in accordance with the fact that these are generally destructive. 

Six-quark decays of W + W~ pairs are detectable but not numerous at E cm = 172 GeV for 
2/cut ~ 0.004, assuming 500 inverse picobarn of luminosity. At the minimum y cnt considered here 
and for the mentioned figure of J Cdt some 8 events should be expected. In contrast, at a LC 
running at the top-antitop threshold, i.e., E cm « 2m t ~ 350 GeV, assuming 100 to 500 fb _1 
per year (e.g., in the TESLA design), one gets between 100 and 500 events, for y cnt = 0.001. 
Presumably, a QCD X-factor of order 1.5-2 should apply to the total production rates of 



e + e~ — > W + W~ — > qqq'q'q"q", in line with the results obtained for e + e~ — > 7*, Z — >• qqq'q' [[36 
so that the actual number of events detected should accordingly be larger. 

As for effects of W + W~ — > 6 quark events onto the line-shape of the mass resonance, 
we have found these negligible. We have performed MINUIT J57] fits of the form 



2 2 

f( m ) = c * 7 2 Xl^ 2 2 + 9(m), (29) 
(m 2 - c 2 ,) 2 + 44 



where the term g(m) is meant to simulate a smooth background due to mis-assigned jets induced 
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by the jet-clustering algorithm, 



g{m) = < 



0. 

c 4 + c 5 (m - c 2 ) + c 6 (m - c 2 ) 2 , (30) 



C4 



l+exp((m-c 5 )/c 6 ) ' 



that is, a null, a three-term polynomial and a smeared step function, see Ref. 0. Clearly, in 
eq. we have assumed a Breit-Wigner shape characterised by a peak height ci, a position c 2 
and a width C3, these corresponding to the normalisation, Mw and Tw, respectively. By adopting 
various algorithms and resolutions, we have never found a difference larger than 10 MeV between 
the c 2 coefficients obtained from various pairs of mass spectra m (one computed with and the 
other without the mentioned correlations), neither at LEP2 nor at a LC. Fig. ^ illustrates typical 
differences, e.g., in the case of an 'average mass'. This can be obtained by applying the Durham 
algorithm with y cnt = 0.001 to the 2 — > 6 process e + e~ — > W + W~ — > qqq'q'q"q" and then forcing 
the six-body final state into a four-body one, by clustering the three softest particles i, j, k into 
one pseudo-particle I with four- momentum pf = pf + Pj + p 1 ^. (the so-called 'E recombination 
scheme' ||). Having done this, one looks at the three possible pairs of di-jet combinations that 
can be formed out of the four surviving four-momenta, rejects the one in which the two most 
energetic particles are put together and plots the average of the other two, M ave . 

We will now proceed to studying the relevance of e + e~ — > W + W~ — > qqq'q'q"q" events in 
the search for new particles at LEP2 and a LC. In doing so, we will use the full ME, with all 
mentioned correlations included. For a start, we have found it rather easy to individuate four 
particle momenta out of the original six in the final state that preserve the typical differences 
between (on the one hand) gluons and (on the other hand) quarks and gluinos in the variables 
which are used to fit the QCD colour factors (see Footnote 1). This can be verified by referring 
to Fig. |], where the differential distributions (normalised to unity) in the four variables^ 

1. the Bengtsson-Zerwas angle xbz] 

2. the (modified) Korner-Schierholz-Willrodt angle $xswi 

3. the (modified) Nachtmann-Reiter angle 9^ R ; 

4. the angle between the two least energetic jets #34; 

have been plotted for e + e~ — > W + W~ — > qqq'q'gg and e + e~ — > W + W~ — > qqq'q'gg events, e.g., 
at LEP2f]. Here, for reference, the gluino mass has been set equal to one (usual jet-finder and 
resolution have been used). The above four-particle quantities have been built by simply using 
the four three-momenta that survive after having removed the pair which yields the invariant 
mass closer to Mw- 

As for the total rates of the various six-jet contributions, these can be found in Fig. |^. The 
CM energies are the same as in Fig. |2|. It is clear from the LEP2 plot that at such a collider 

8 See, e.g., Ref. [[38| for the definition of the angles and for some typical spectra in the case of e + e~ — > 4 parton 
processes at the Z peak. 
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Those for e+e" -> W+W~ -> qqq'q'q"q" 
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there is no chance of selecting a statistically significant sample of SUSY events. If one assumes, 
say, 500 inverse picobarn to be collected at LEP2, then, from looking at the upper part of Fig. || 
it follows that the number of SUSY events produced should be between 1 and 8, for y cut = 0.001 
in the Durham algorithm, depending on the actual value of the gluino mass ! In contrast, 
one should expect a LC to be an excellent laboratory for gluino searches in W + W~ decays. 
(Incidentally, notice that at a LC the cross section for W + W~ production is even larger than 
that for the Z.) In fact, although the LC production rate of e + e~ — > W + W~ is smaller than 
the LEP2 one, the instantaneous luminosity is in contrast much higher, since some 100 to 500 
fb _1 of data per annum are expected to be collected. If one goes back to the lower part of Fig. |] 
and considers a LC running at 350 GeV, then one should expect to produce approximately 120 
to 600 events per year with very light gluinos, for y cut = 0.001 in the Durham scheme (possibly, 
twice as much, accounting for the i^-factor). 

Varying the CM energy of the colliding e + e~ beams has little effects on the effectiveness of 
our simple procedure based on the M\y selection of four four-momenta. This can be seen by 
defining the Lorentz-invariant (contrary to angles) quantity 

v = Nv e(n,P2,P 3 ,P,) 
s 2 

where the numerator represents the contraction of the jet four-momenta Pi,P2,P3,Pl used in 
Fig. | with the Levi-Civita tensor e^ upi j of Sect. |3] and where the factor Ny has been introduced 
for scaling purposes. The quantity in eq. (|3T| ) is presented in Fig. || e.g., for E cm = 172 GeV, 
using Ny = 1000: compare to Fig. |], where the same energy was used. In this case, we again 
have looked at the four-quark-two-gluon and four-quark-two-gluino final states only, as the six- 
quark contribution behaves rather similarly to the SUSY one. Basically, V is a measure of the 
acoplanarity of the event this in turn quantifying the relative orientation of the planes 



spanned by, on the one hand, the two most energetic particles and, on the other hand, the two 
least energetic ones. In fact, it should be recalled that the first three angles introduced above are 
nothing else than a different way of describing the helicity property that in a g* — > gg splitting 
the two gluons (they are spin 1 bosons) tend to lie in the same plane of the two quarks which 
originally emitted the virtual gluon, whereas in g* — > qq and g* — > gg splittings the two quarks 
and the two gluinos (they are both spin 1/2 fermions) tend to be in a perpendicular one. Fig. |6] 
eloquently confirms this dynamics. 

Under these circumstances then, dedicated analyses in the angular variables of Fig. f| could 
well be attempted at a LC (also notice the somewhat improved 'SUSY to ordinary QCD' pro- 
duction ratio respect to LEP2). To select a six-jet sample from W + W~ decays should be rather 
straightforward, we believe, by removing those where none of the di-jet invariant masses repro- 
duce M w , most of which would come from 0{af) QCD events p5] , and those in which one or 
more three-jet masses reconstruct m t = 175 GeV, as it occurs in single- and double-top events 
POfl . In principle, one could also resort to semileptonic decays, i.e., W + W~ — > 2 jets £ ± plus 
missing energy. In practise, though, the loss of kinematic constraints, because of the neutrino 
escaping detection, would render the background suppression less effective. Furthermore, to 
exploit the sort of vertex tagging procedure described in Ref. |2^j and recalled in Sect. ^| could 



be even more fruitful, given the high efficiency and purity foreseen for such a technique by 
the time the electron-positron linear colliders will have started operation, but provided that 
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t 5 ~ 2 — 5- ~ Tb (e.g., in the SUGRA scenario), with t& ~ 1.6 • 10 12 sec and where M 5 is the 

typical squark mass. 

We now turn our attention to the case of Higgs searches in the 66 + 4 jet channel, at a LC0. 
Though e + e~ — > W + W~ — > qqq'q'bb events can be relevant as background processes in both 
the SM and the MSSM, for reasons of space, we illustrate here the phenomenology of the latter 
model only. The former case can easily be dealt with by the reader itself, by referring to the 
specialised bibliography on the subject [[41]]. 

Fig. ^presents the Higgs rates at a LC with E cm = 350 (above) and 500 (below) GeV, for two 
reference values of tan/?, 3.0 and 30., as a function of the scalar Higgs masses. The background 
considered here clearly does not depend on either of them, so it is simply indicated by an arrow 
in both plots, in correspondence of the transition value of the Higgs mass between the light, M h , 
and heavy, Mh, regime. As usual, six-jet final states (here, with two 6 quarks) are selected using 
the Durham jet-finder with cut-off y cut = 0.001. For such jet selection enforced, the background 
overwhelms both scalar Higgs signals, e + e~ — > <3?Z and e + e~ — > $v4, over a large interval in 
and Mh- 

However, the situation is in reality less dramatic than it would appear from Fig. [7] only, if 
one refers to Fig. |] too. In fact, whereas all the decay products emerging from the W + W~Z 
and W + W~A intermediate states of the Higgs signals are naturally energetic and far apart, 
the 6 quarks generated in W + W~ events tend to be soft and collinear, owning to the dominant 
infrared dynamics of the gluon splitting (see the dashed lines in the top plots of Figs. [TO] and 
|I~2| below). In other terms, whereas to increase the value of y cnt would affect both signals only 
slightly, this is no longer true for the background. In fact, in one increases the resolution, 
e.g., by a factor of five, to y cut = 0.005, the latter decreases by a factor of about 90(25) at 
E C m = 350(500) GeV. In correspondence, the typical loss for each of the former is less than a 
factor of 5 at both energies. 

Therefore, even in those cases where the background is apparently well above the signal, 
e.g, for M A = 220 GeV and tan/3 = 3.0 GeV (corresponding to M H ~ 230 GeV), the latter 
being dominated at 350 GeV by HZ — > W + W~ Z production and decay (the asterisk in the 
upper frame of Fig. |7|), a judicious choice of y cnt combined with a dedicate selection in Mij (with 
i < j = 1, ...4, the six di-jet invariant masses that can be reconstructed from the light quark jets 
ordered in energy, see Fig. |9|) and M b i (that of the two 6 jets, see top frame of Fig. [TUp around 
the W and Z masses, respectively, might allow for the remotion of e + e~ — > W + W~ — > qqq'q'bb 
events in the spectrum of the four-light- quark invariant mass, in which the Higgs peak generated 
in the production and decay sequence e + e~ — > HZ — > W + W~Z — > qqq'q'bb should appear 
(bottom frame of Fig. [ll]). (Notice that bins are there two GeV wide, being the expected 
M 4 g resolution more realistically, say, five times as large, so that the signal would actually be 
well below the background in this distribution, after jet-selection cuts only.) However, even for 
/ Cdt = 500 fb _1 , the event rate is rather poor in this case, about four events per year. 

In contrast, for other settings of the MSSM parameter spaces, the signal would clearly be 

10 Notice that in all forthcoming plots we have not included a multiplicative factor e^, accounting for the finite 
efficiency of tagging the two heavy quarks. We assume th to be large enough so that the reducible background 
from e + e~ — > W + W~ — > qqq'q'gg studied in Ref. JToJ and in the first part of this Section can easily be filtered 
out of the bb + 4 jet sample. 
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visible above the QCD noise considered here, even before the implementation of the Higgs 
selection cuts. As illustrative example, we consider again the same point, Ma = 220 GeV 
and tan/3 = 3.0 GeV (Mjj « 230 GeV), but at 500 GeV, where the dominant Higgs channel 
now involves HA — > W + W~A production and decay (see the asterisk in the bottom plot of 
Fig. [?[). For such a choice, the My/ resonance is already much higher for the signal than for the 
background (see Fig. 11], where the arrows denote the height of the peaks for the latter) and 
both the A and H Breit-Wigner shapes clearly stick out in the M b i and M 4(? spectra, respectively 
(see Fig. [12]), even if the four-jet mass resolution is much larger than 2 GeV. In this case, the 
yearly production rate of the signal would be 140 events (again, assuming / Cdt = 500 fb _1 ). 

The discussion for the case of the light scalar Higgs of the MSSM is rather similar, so we do 
not repeat it here. In this case, in general, the MSSM parameter space accessible via the Higgs 
signature h —>■ W + W~ — ► 4 jets is much reduced though, as only a light scalar with mass at the 
very upper hand of its allowed range can decay in such a channel (in which one of the two 
bosons is off-shell): see Fig. [7|. 



5 Conclusions 

In the end, 0(a 2 s ) decays into six-fermions of W + W~ pairs produced in e + e~ scatterings can 
be detected at LEP2 and a LC as well, the latter with 350 GeV S E cm S 500 GeV. At both 
colliders though, they have little relevance in My/ measurements. As for New Physics analyses, 
such events can be important in a LC environment. Firstly, they can produce very light gluinos 
at statistically significant rate. These sparticles have survived the LEP constraints and in some 
SUSY scenarios could well be the next-to-lightest ones, all other being much heavier. The new 
SUSY signals could be searched for in six-jet samples in which only two of the jets reconstruct the 
W mass. Using the remaining four-jet subset one could either fit the QCD colour factors to the 
shape of some typical angular distributions or exploit the tagging of a displaced vertex. Secondly, 
they can represent an overwhelming background in the search for Higgs bosons, both in the SM 
and in the MSSM, produced via e+e" -> Higgs Z -> W + W~Z and e+e~ -> <M -> W + W~A, 
where Higgs = 0, $ and $ = H,h are the scalar particles and in which W + W~ — > 4 (light) jets. 
However, to cut around the Z and/or A masses in the bb subsystem, should in general allow 
one to reduce drastically such a QCD noise in the Higgs candidate sample of six-jets with two 
displaced vertices, provided a high b tagging efficiency can be achieved. 

The calculation of the exact matrix element for the 2 — > 6 process e + e~ — ► W + W~ — > 
qqq'q'q"q" through the order 0{a^ m a 2 s ) has been performed by means of helicity amplitude 
methods. This has led to compact analytic expressions of the MEs and to their fast implemen- 
tation in a FORTRAN code which can profitably be exploited in high statistic MC simulations. 
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Diagrams by MadGraph 




graph 7 graph 8 



Figure 1: Feynman diagrams contributing to e + (l) + e _ (2) — > g(3) + g(4) + g'(5) + <f (6) + 
q"(7) + <?" (8), via W + W~ production and decay, through the order 0{a^ m a 2 s ). Here, we assume 
q" 7^ q', q. An internal wavy line represents a 7, Z or a W ± , as appropriate. 
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Figure 2: Cross section for e + e~ — > W + W~~ — > qqq'q'q"q" events as a function of the resolution 
parameter y cut in the Durham jet-finder, at E cm = 172 (above) and 350 (below) GeV. The 
summation over all possible combinations of flavours q, q' and q" has been performed. Solid: 
without Fermi-Dirac correlations. Dashed: with Fermi-Dirac correlations. 



e e -> W W -> 6q 



0.004 



0.003 



0.002 



0.001 



0.000 



0.25 



0.20 



0.15 



0.10 



0.05 



0.00 



20 




40 60 

M ave (GeV) 

+ - + - 
e e -> W W -» 6q 




(GeV) 



Figure 3: Differential distributions in the 'average value' of the best reconstructed W ± masses 
for e + e~ — > W + W~ — ► qqq'q'q"q" events (solid and dashed lines) and a for pure phase-space 
model of six-particle production (dotted line, below only), with y cut = 0.001 in the Durham 
jet-finder, at E cm = 172 GeV. The summation over all possible combinations of flavours q, q' 
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Figure 4: Differential distributions in the cosine of the following angular variables: (top-left) xbz, 
(top-right) $ksw> (bottom-left) 6^ R and (bottom-right) 6> 34 , for e + e~ — > W + W~ — > qqq'q'gg 
(solid lines) and e + e~ — > W + W~ — > qqq'q'gg (dashed lines) events, for = 1 GeV, with 
2/cut = 0.001 in the Durham jet-finder, at E cm = 172 GeV. The summation over all possible 
combinations of flavours q and q' has been performed. Spectra are normalised to unity. 
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Figure 5: Cross section for e + e~ — > W + W~ — > qqq'q'gg (solid lines), e + e~ — > W + W~ — > 
qqq'q'q"q" (dashed lines) and e + e~ — > W + W~ — > qqq'q'gg (dotted lines) events as a function of 
the gluino mass, with y cut = 0.001 in the Durham jet-finder, at _E cm = 172 (upper plot) and 350 
GeV (lower plot). The summation over all possible combinations of flavours q and q' has been 
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Figure 6: Differential distributions in the quantity defined in eq. fl3lD of the resonant four-jet 
subsystem for e + e~ — ► W + W~ — > qqq'q'gg (solid lines) and e + e~ — > W + W~~ — > qqq'q'gg (dashed 
lines) events, for rrig = 1 GeV, with y cut = 0.001 in the Durham jet-finder, at E cm = 172 GeV. 
The summation over all possible combinations of flavours q and q' has been performed. Spectra 
are normalised to unity. 
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Figure 7: Cross section for e+e~ — > $Z — > W + W~Z — > qqq'q'bb, e + e~ — >■ <M — > W + W~A — > 
qqq'q'bb (solid, dashed, dotted and dot-dashed lines) and e + e~~ — > W + W~ — > qqq'q'bb (arrows) 
events as a function of the scalar Higgs masses M h and M#, with y cut = 0.001 in the Durham 
jet-finder, at _E cm = 350 (upper plot) and 500 (lower plot) GeV. The summation over all possible 
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Figure 8: Cross section for e + e~ — > W + W~ — > qqq'q'bb events as a function of the resolution 
parameter y cut in the Durham jet-finder, at E cm = 350 (solid line) and 500 (dashed line) GeV. 
The summation over all possible combinations of flavours q and q' has been performed. 
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Figure 9: Differential distributions in the invariant mass of the energy-ordered di-jet pairs, 
Mij with i < j = 1, ...4, for e + e~ — > iJZ — > W + W~Z — > qqq'q'bb (solid lines) and e + e~ — > 
W + W~ _,. qqq'q'bb (dashed lines) events, in the MSSM for M A = 220 GeV and tan/5 = 3.0, 
with y cnt = 0.001 in the Durham jet-finder, at E cm = 350 GeV. The summation over all possible 
combinations of flavours q and q' has been performed. Spectra are normalised to the total cross 
sections. 
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Figure 10: Differential distributions in the invariant mass of the bb, M b i (upper plot), and four- 
light-quark, M iq (lower plot), systems, for e + e~ — > HZ — > W + W~Z — > qqq'q'bb (solid lines) 
and e + e~ -> W+Vy - -> (dashed lines) events, in the MSSM for M A = 220 GeV and 

tan/3 = 3.0, with y cut = 0.001 in the Durham jet-finder, at E cm = 350 GeV. The summation 
over all possible combinations of flavours q and q' has been performed. Spectra are normalised 
to the total cross sections. 
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Figure 11: Differential distributions in the invariant mass of the energy-ordered di-jet pairs, 
Mij with i < j = 1, ...4, for e + e~ — > iJA — > W + W~A — > qqq'q'bb (solid lines) and e + e~ — > 
W + W~ — >■ qqq'q'bb (dashed lines, where the peak height is denoted by an arrow) events, in the 
MSSM for M A = 220 GeV and tan (5 = 3.0, with y cut = 0.001 in the Durham jet-finder, at 
E cm = 500 GeV. The summation over all possible combinations of flavours q and q' has been 
performed. Spectra are normalised to the total cross sections. 
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Figure 12: Differential distributions in the invariant mass of the bb, M b i (upper plot), and four- 
light-quark, M iq (lower plot), systems, for e + e~ — > HA — > W + W~A — > qqq'q'bb (solid lines) 
and e + e~ -> W + W~ -> qqq'q'bb (dashed lines) events, in the MSSM for M A = 220 GeV and 
tan/3 = 3.0, with y cut = 0.001 in the Durham jet-finder, at E cm = 500 GeV. The summation 
over all possible combinations of flavours q and q' has been performed. Spectra are normalised 
to the total cross sections. 



